Regulation of long-range BMP gradients and embryonic polarity by propagation of local calcium-firing activity.

Many amniote vertebrate species including humans can form identical twins from a single embryo, but this only occurs rarely. It has been suggested that the primitive-streak-forming embryonic region emits signals that inhibit streak formation elsewhere but the signals involved, how they are transmitted and how they act has not been elucidated. Here we show that short tracks of calcium firing activity propagate through extraembryonic tissue via gap junctions and prevent ectopic primitive streak formation in chick embryos. Cross-regulation of calcium activity and an inhibitor of primitive streak formation (Bone Morphogenetic Protein, BMP) via NF-κB and NFAT establishes a long-range BMP gradient spanning the embryo. This mechanism explains how embryos of widely different sizes can maintain positional information that determines embryo polarity. We provide evidence for similar mechanisms in two different human embryo models and in Drosophila, suggesting an ancient evolutionary origin.

A gene regulatory network for neural induction.

During early vertebrate development, signals from a special region of the embryo, the organizer, can redirect the fate of non-neural ectoderm cells to form a complete, patterned nervous system. This is called neural induction and has generally been imagined as a single signalling event, causing a switch of fate. Here, we undertake a comprehensive analysis, in very fine time course, of the events following exposure of competent ectoderm of the chick to the organizer (the tip of the primitive streak, Hensen's node). Using transcriptomics and epigenomics we generate a gene regulatory network comprising 175 transcriptional regulators and 5614 predicted interactions between them, with fine temporal dynamics from initial exposure to the signals to expression of mature neural plate markers. Using in situ hybridization, single-cell RNA-sequencing, and reporter assays, we show that the gene regulatory hierarchy of responses to a grafted organizer closely resembles the events of normal neural plate development. The study is accompanied by an extensive resource, including information about conservation of the predicted enhancers in other vertebrates.

The embryonic node behaves as an instructive stem cell niche for axial elongation.

In warm-blooded vertebrate embryos (mammals and birds), the axial tissues of the body form from a growth zone at the tail end, Hensen's node, which generates neural, mesodermal, and endodermal structures along the midline. While most cells only pass through this region, the node has been suggested to contain a small population of resident stem cells. However, it is unknown whether the rest of the node constitutes an instructive niche that specifies this self-renewal behavior. Here, we use heterotopic transplantation of groups and single cells and show that cells not destined to enter the node can become resident and self-renew. Long-term resident cells are restricted to the posterior part of the node and single-cell RNA-sequencing reveals that the majority of these resident cells preferentially express G2/M phase cell-cycle-related genes. These results provide strong evidence that the node functions as a niche to maintain self-renewal of axial progenitors.

Molecular anatomy of the pre-primitive-streak chick embryo.

The early stages of development of the chick embryo, leading to primitive streak formation (the start of gastrulation), have received renewed attention recently, especially for studies of the mechanisms of large-scale cell movements and those that position the primitive streak in the radial blastodisc. Over the long history of chick embryology, the terminology used to define different regions has been changing, making it difficult to relate studies to each other. To resolve this objectively requires precise definitions of the regions based on anatomical and functional criteria, along with a systematic molecular map that can be compared directly to the functional anatomy. Here, we undertake these tasks. We describe the characteristic cell morphologies (using scanning electron microscopy and immunocytochemistry for cell polarity markers) in different regions and at successive stages. RNAseq was performed for 12 regions of the blastodisc, from which a set of putative regional markers was selected. These were studied in detail by in situ hybridization. Together this provides a comprehensive resource allowing the community to define the regions unambiguously and objectively. In addition to helping with future experimental design and interpretation, this resource will also be useful for evolutionary comparisons between different vertebrate species.

A 3D molecular atlas of the chick embryonic heart.

We present a detailed analysis of gene expression in the 2-day (HH12) embryonic chick heart. RNA-seq of 13 micro-dissected regions reveals regionalised expression of 15,570 genes. Of these, 132 were studied by in situ hybridisation and a subset (38 genes) was mapped by Optical Projection Tomography or serial sectioning to build a detailed 3-dimensional atlas of expression. We display this with a novel interactive 3-D viewer and as stacks of sections, revealing the boundaries of expression domains and regions of overlap. Analysis of the expression domains also defines some sub-regions distinct from those normally recognised by anatomical criteria at this stage of development, such as a previously undescribed subdivision of the atria into two orthogonal sets of domains (dorsoventral and left-right). We also include a detailed comparison of expression in the chick with the mouse and other species.

BRepertoire: a user-friendly web server for analysing antibody repertoire data.

Antibody repertoire analysis by high throughput sequencing is now widely used, but a persisting challenge is enabling immunologists to explore their data to discover discriminating repertoire features for their own particular investigations. Computational methods are necessary for large-scale evaluation of antibody properties. We have developed BRepertoire, a suite of user-friendly web-based software tools for large-scale statistical analyses of repertoire data. The software is able to use data preprocessed by IMGT, and performs statistical and comparative analyses with versatile plotting options. BRepertoire has been designed to operate in various modes, for example analysing sequence-specific V(D)J gene usage, discerning physico-chemical properties of the CDR regions and clustering of clonotypes. Those analyses are performed on the fly by a number of R packages and are deployed by a shiny web platform. The user can download the analysed data in different table formats and save the generated plots as image files ready for publication. We believe BRepertoire to be a versatile analytical tool that complements experimental studies of immune repertoires. To illustrate the server's functionality, we show use cases including differential gene usage in a vaccination dataset and analysis of CDR3H properties in old and young individuals. The server is accessible under http://mabra.biomed.kcl.ac.uk/BRepertoire.

Promiscuous antibodies characterised by their physico-chemical properties: From sequence to structure and back.

Human B cells produce antibodies, which bind to their cognate antigen based on distinct molecular properties of the antibody CDR loop. We have analysed a set of 10 antibodies showing a clear difference in their binding properties to a panel of antigens, resulting in two subsets of antibodies with a distinct binding phenotype. We call the observed binding multiplicity 'promiscuous' and selected physico-chemical CDRH3 characteristics and conformational preferences may characterise these promiscuous antibodies. To classify CDRH3 physico-chemical properties playing a role in their binding properties, we used statistical analyses of the sequences annotated by Kidera factors. To characterise structure-function requirements for antigen binding multiplicity we employed Molecular Modelling and Monte Carlo based coarse-grained simulations. The ability to predict the molecular causes of promiscuous, multi-binding behaviour would greatly improve the efficiency of the therapeutic antibody discovery process.

PinSnps: structural and functional analysis of SNPs in the context of protein interaction networks.

UNLABELLED: We present a practical computational pipeline to readily perform data analyses of protein-protein interaction networks by using genetic and functional information mapped onto protein structures. We provide a 3D representation of the available protein structure and its regions (surface, interface, core and disordered) for the selected genetic variants and/or SNPs, and a prediction of the mutants' impact on the protein as measured by a range of methods. We have mapped in total 2587 genetic disorder-related SNPs from OMIM, 587 873 cancer-related variants from COSMIC, and 1 484 045 SNPs from dbSNP. All result data can be downloaded by the user together with an R-script to compute the enrichment of SNPs/variants in selected structural regions. AVAILABILITY AND IMPLEMENTATION: PinSnps is available as open-access service at http://fraternalilab.kcl.ac.uk/PinSnps/ CONTACT: franca.fraternali@kcl.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Anatomy of protein disorder, flexibility and disease-related mutations.

Integration of protein structural information with human genetic variation and pathogenic mutations is essential to understand molecular mechanisms associated with the effects of polymorphisms on protein interactions and cellular processes. We investigate occurrences of non-synonymous SNPs in ordered and disordered protein regions by systematic mapping of common variants and disease-related SNPs onto these regions. We show that common variants accumulate in disordered regions; conversely pathogenic variants are significantly depleted in disordered regions. These different occurrences of pathogenic and common SNPs can be attributed to a negative selection on random mutations in structurally highly constrained regions. New approaches in the study of quantitative effects of pathogenic-related mutations should effectively account for all the possible contexts and relative functional constraints in which the sequence variation occurs.

AP1S3 mutations are associated with pustular psoriasis and impaired Toll-like receptor 3 trafficking.

Adaptor protein complex 1 (AP-1) is an evolutionary conserved heterotetramer that promotes vesicular trafficking between the trans-Golgi network and the endosomes. The knockout of most murine AP-1 complex subunits is embryonically lethal, so the identification of human disease-associated alleles has the unique potential to deliver insights into gene function. Here, we report two founder mutations (c.11T>G [p.Phe4Cys] and c.97C>T [p.Arg33Trp]) in AP1S3, the gene encoding AP-1 complex subunit σ1C, in 15 unrelated individuals with a severe autoinflammatory skin disorder known as pustular psoriasis. Because the variants are predicted to destabilize the 3D structure of the AP-1 complex, we generated AP1S3-knockdown cell lines to investigate the consequences of AP-1 deficiency in skin keratinocytes. We found that AP1S3 silencing disrupted the endosomal translocation of the innate pattern-recognition receptor TLR-3 (Toll-like receptor 3) and resulted in a marked inhibition of downstream signaling. These findings identify pustular psoriasis as an autoinflammatory phenotype caused by defects in vesicular trafficking and demonstrate a requirement of AP-1 for Toll-like receptor homeostasis.

The Protein 4.1 family: hub proteins in animals for organizing membrane proteins.

Proteins of the 4.1 family are characteristic of eumetazoan organisms. Invertebrates contain single 4.1 genes and the Drosophila model suggests that 4.1 is essential for animal life. Vertebrates have four paralogues, known as 4.1R, 4.1N, 4.1G and 4.1B, which are additionally duplicated in the ray-finned fish. Protein 4.1R was the first to be discovered: it is a major mammalian erythrocyte cytoskeletal protein, essential to the mechanochemical properties of red cell membranes because it promotes the interaction between spectrin and actin in the membrane cytoskeleton. 4.1R also binds certain phospholipids and is required for the stable cell surface accumulation of a number of erythrocyte transmembrane proteins that span multiple functional classes; these include cell adhesion molecules, transporters and a chemokine receptor. The vertebrate 4.1 proteins are expressed in most tissues, and they are required for the correct cell surface accumulation of a very wide variety of membrane proteins including G-Protein coupled receptors, voltage-gated and ligand-gated channels, as well as the classes identified in erythrocytes. Indeed, such large numbers of protein interactions have been mapped for mammalian 4.1 proteins, most especially 4.1R, that it appears that they can act as hubs for membrane protein organization. The range of critical interactions of 4.1 proteins is reflected in disease relationships that include hereditary anaemias, tumour suppression, control of heartbeat and nervous system function. The 4.1 proteins are defined by their domain structure: apart from the spectrin/actin-binding domain they have FERM and FERM-adjacent domains and a unique C-terminal domain. Both the FERM and C-terminal domains can bind transmembrane proteins, thus they have the potential to be cross-linkers for membrane proteins. The activity of the FERM domain is subject to multiple modes of regulation via binding of regulatory ligands, phosphorylation of the FERM associated domain and differential mRNA splicing. Finally, the spectrum of interactions of the 4.1 proteins overlaps with that of another membrane-cytoskeleton linker, ankyrin. Both ankyrin and 4.1 link to the actin cytoskeleton via spectrin, and we hypothesize that differential regulation of 4.1 proteins and ankyrins allows highly selective control of cell surface protein accumulation and, hence, function. This article is part of a Special Issue entitled: Reciprocal influences between cell cytoskeleton and membrane channels, receptors and transporters. Guest Editor: Jean Claude Hervé

Mapping disease-related missense mutations in the immunoglobulin-like fold domain of lamin A/C reveals novel genotype-phenotype associations for laminopathies.

Mutations in A-type nuclear lamins cause laminopathies. However, genotype-phenotype correlations using the 340 missense mutations within the LMNA gene are unclear: partially due to the limited availability of three-dimensional structure. The immunoglobulin (Ig)-like fold domain has been solved, and using bioinformatics tools (including Polyphen-2, Fold X, Parameter OPtimized Surfaces, and PocketPicker) we characterized 56 missense mutations for position, surface exposure, change in charge and effect on Ig-like fold stability. We find that 21 of the 27 mutations associated with a skeletal muscle phenotype are distributed throughout the Ig-like fold, are nonsurface exposed and predicted to disrupt overall stability of the Ig-like fold domain. Intriguingly, the remaining 6 mutations clustered, had higher surface exposure, and did not affect stability. The majority of 9 lipodystrophy or 10 premature aging syndrome mutations also did not disrupt Ig-like fold domain stability and were surface exposed and clustered in distinct regions that overlap predicted binding pockets. Although buried, the 10 cardiac mutations had no other consistent properties. Finally, most lipodystrophy and premature aging mutations resulted in a -1 net charge change, whereas skeletal muscle mutations caused no consistent net charge changes. Since premature aging, lipodystrophy and the subset of 6 skeletal muscle mutations cluster tightly in distinct, charged regions, they likely affect lamin A/C -protein/DNA/RNA interactions: providing a consistent genotype-phenotype relationship for mutations in this domain. Thus, this subgroup of skeletal muscle laminopathies that we term the 'Skeletal muscle cluster', may have a distinct pathological mechanism. These novel associations refine the ability to predict clinical features caused by certain LMNA missense mutations.

Protein-protein interaction networks studies and importance of 3D structure knowledge.

Protein-protein interaction networks (PPINs) are a powerful tool to study biological processes in living cells. In this review, we present the progress of PPIN studies from abstract to more detailed representations. We will focus on 3D interactome networks, which offer detailed information at the atomic level. This information can be exploited in understanding not only the underlying cellular mechanisms, but also how human variants and disease-causing mutations affect protein functions and complexes' stability. Recent studies have used structural information on PPINs to also understand the molecular mechanisms of binding partner selection. We will address the challenges in generating 3D PPINs due to the restricted number of solved protein structures. Finally, some of the current use of 3D PPINs will be discussed, highlighting their contribution to the studies in genotype-phenotype relationships and in the optimization of targeted studies to design novel chemical compounds for medical treatments.

Specialized Dynamical Properties of Promiscuous Residues Revealed by Simulated Conformational Ensembles.

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein-protein interaction prediction and design methods.

B-RAF mutant alleles associated with Langerhans cell histiocytosis, a granulomatous pediatric disease.

BACKGROUND: Langerhans cell histiocytosis (LCH) features inflammatory granuloma characterised by the presence of CD1a+ dendritic cells or 'LCH cells'. Badalian-Very et al. recently reported the presence of a canonical (V600E)B-RAF mutation in 57% of paraffin-embedded biopsies from LCH granuloma. Here we confirm their findings and report the identification of two novel B-RAF mutations detected in LCH patients. METHODS AND RESULTS: Mutations of B-RAF were observed in granuloma samples from 11 out of 16 patients using 'next generation' pyrosequencing. In 9 cases the mutation identified was (V600E)B-RAF. In 2 cases novel polymorphisms were identified. A somatic (600DLAT)B-RAF insertion mimicked the structural and functional consequences of the (V600E)B-RAF mutant. It destabilized the inactive conformation of the B-RAF kinase and resulted in increased ERK activation in 293 T cells. The (600DLAT)B-RAF and (V600E)B-RAF mutations were found enriched in DNA and mRNA from the CD1a+ fraction of granuloma. They were absent from the blood and monocytes of 58 LCH patients, with a lower threshold of sequencing sensitivity of 1%-2% relative mutation abundance. A novel germ line (T599A)B-RAF mutant allele was detected in one patient, at a relative mutation abundance close to 50% in the LCH granuloma, blood monocytes and lymphocytes. However, (T599A)B-RAF did not destabilize the inactive conformation of the B-RAF kinase, and did not induce increased ERK phosphorylation or C-RAF transactivation. CONCLUSIONS: Our data confirmed presence of the (V600E)B-RAF mutation in LCH granuloma of some patients, and identify two novel B-RAF mutations. They indicate that (V600E)B-RAF and (600DLAT)B-RAF mutations are somatic mutants enriched in LCH CD1a(+) cells and absent from the patient blood. Further studies are needed to assess the functional consequences of the germ-line (T599A)B-RAF allele.